Spin interactions of interstitial Mn ions in ferromagnetic GaMnAs 
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The recMitly reported Rutherford backscattering and particle-induced X-ray emission 
experimentaj have revealed that in low-temperature MBE grown Gai-a^Mn^^As a significant part 
of the incorporated Mn atoms occupies tetrahedral interstitial sites in the lattice. Here we study 
the magnetic properties of these interstitial (Mni) ions. We show that they do not participate in 
the hole-induced ferromagnetism. Moreover, Mni double donors may form pairs with the nearest 
substitutional (MuQa) acceptors - our calculations evidence that the spins in such pairs are anti- 
ferromagnetically coupled by the superexchange. We also show that for the Mn ion in the other, 
hexagonal, interstitial position (which seems to be the case in the Gai-x-yMnxBeyAs samples) the 
p-d interactions with the holes, responsible for the ferromagnetism, are very much suppressed. 



The incorporation of transition metal ions into the 
III-V host semiconductors by low-temperature molecu- 
lar beam epitaxy (LT MBE), i.e., the discovery of ferro- 
magnetic dilute magnetic semiconductors (DMS) in the 
pioneering work by Munekata et alB, was a major step 
towards the integration of the spin degrees of freedom 
with the semiconducting properties in the same mate- 
rial. Still, the prospects for practical applications of DMS 
in "spintronic" devices depend crucially on the possibil- 
ities to increase in these materials the temperature of 
the transition to the ferromagnetic phase. The highest 
Curie temperatures (Tc) in DMS have been obtained by 
a substitution of Mn for Ga in GaAs, which was com- 
plemented by a post-growth annealing in temperatures 
only slightly exceeding the LT MBE growth tempera- 
ture. Until recently the Tcpf4rl|Q K seemed to be the 
upper limit for this materialElElffl. In the last months, 
however, considerably higher values of Tc in GaMnAs, 
even exceeding 15fl K for , thin films, have been reported 
by several groupsall3llilll3. This progress has been made 
basically by an optimization of the annealing time and 
temperature. 

In the theoretical models describing the ferromag- 
netism in DMS (e.g., in Ref. |l|,^jl5|) the Tc is expected 
to increase with both, the magnetic ions and hole con- 
centrations. In LT MBE grown Gai-xMn^As this was 
indeed experimentally-established for Mn concentrations 
up to about x=0. 07,13. The Mn ion in the substitu- 
tional position in the GaAs lattice (Mnca) acts as an 
acceptor, but in all Gai_xMnxAs samples the hole con- 
centration is substantially lower than the Mn content. 
This has been ascribed to the presence of compensating 
donors, in particular to the formation of arsenic anti- 
sites (Asca) durins-the epitaxial growth of Gai_xMnxAs 
at As overpressuraiBEZl. In Refj^,^ and |l^ the observed 
annealing-induced changes of the Tc were attributed 
solely to the decrease of the concentration of arsenic 
antisites leading to the increase of the hole concentra- 
tion. These antisites, however, are relatively stable de- 
fects - it was shown that to remove Asca from LT MBE 



grown GaAs the annealing temperatures above 450 C 
are neededtEl. Recently, simultaneous channeling Ruther- 
ford backscattering (c-RBS) and particle-induced X-ray 
emission-. (c-PIXE) experiments shed new light on this 
problcmtl Namely, they have revealed that in LT MBE 
grown ferromagnetic Gai.xMuxAs with high x a signif- 
icant fraction of incorporated Mn atoms (ca 15% for 
the as-grown Gao.91Mno.09As sample) occupies well de- 
fined, commensurate with the GaAs lattice interstitial 
positions. 

In the diamond cubic crystal lattice there are two pos- 
sible interstitial positions, the so called tetrahedral and 
hexagonal sites, in which the atoms are shadowed along 
(100) and (111) direction and exposed in the (110) axial 
channel, as seen at the experiment. They can be dis- 
tinguished by jStudying angular scans around the (110) 
axial directioncJ. The scans presented in Ref. |l| sug- 
gested that the interstitial Mn ions (Mnj) observed in 
Gai.xMuxAs occupy the tetrahedral sites, in which the 
interstitial is surrounded by four nearest neighbors, as 
presented in Fig. ^ 

The Mni serve, like Asca, as double donors, decreasing 
the hole concentration. The results presented in Ref. |l| 
directly showed that in the process of LT annealing the 
Mni ions are moved to random, incommensurate with the 
GaAs lattice positions (e.g., MnAs clusters), in which 
the Mn ions are electrically inactive. Thus, in the an- 
nealed samples the concentration of the compensating 
Mni donors decreases considerably whereas the hole con- 
centration increases and the observed Curie temperature 
is much higher. Moreover, it was demonstrated that the 
appropriate annealing increases the saturation magneti- 
zation, i.e., that ptlie. jpresence of Mni reduces the net 
magnetic momentcrna'u. 

The described above experimental results stimulated 
theoretical studies on the formation and properties of in- 
terstitial Mn in the GaMnAs ternary compound. First, 
the electronic structure of the GaMnAs with Mn in sub- 
stitutional and interstitial position was calculated by 
ab initio methods, showing that indeed Mn interstitials 
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FIG. 1: The nearest four cation and six anion neighbors for an 
ion in the tetrahedral interstitial position in the zinc-blende 
lattice. 



act as double donorscil. In a recent papeicj the self- 
compensation of Mn in such semiconductors was studied 
within the density-functional theory. In Ref. it was 
shown that interstitial Mn can be easily formed near the 
surface. Here we consider the spin properties of intersti- 
tial magnetic ions. We study the spin interactions for the 
Mni ion in the tetraheadral interstitial position in order 
to provide theoretical basis for the understanding of the 
experimental findings concerning the magnetic behavior 
of the as-grown and annealed Gai_xMnxAs samples. We 
analyze the hybridization of the d-orbitals of these ions 
with the valence band p-states. This effect is essential 
for both, the superexchange and the RKKY-type, dom- 
inant ion-ion interactions in DMS. It is widely accepted 
that the latter mechanism is responsible for the hole- 
induced ferromagnetism in III-V DMS and that the Tc 
depends crucially on the p-d hybridization - within the 
Zener modeluJ Tc is proportional to the square of the 
kinetic p-d exchange constant /3, i.e., to the fourth power 
of the hybridization constant V at the centre of the Bril- 
louin zone. 

The valence band states in Gai.xMuxAs are built pri- 
marily from the anion p-orbitals, thus the p-d hybridiza- 
tion for a given magnetic ion is determined by the posi- 
tions of its nearest-neighbor anions. In zinc-blende lattice 
of GaAs, the Mn ion in the cation substitutional position 
has four anion nearest neighbors at the distances a\/3/4 
(where a is the lattice constant) along the [1, 1, 1], [1,-1,- 
1], [-1, 1,-1] and [-1,-1, 1] directions. For these positions 
the inter-atomic matrix elements, E^^xy, E^^yz, E^^zx, 
etc., exprjessed in terms of the Harrison parameters Vpda- 
and V^(j7rE3, add up constructively to the hybridization 
constant V in the hybridization Hamiltionian , with 
different weights for different points of the Brillouin zone. 
At the point fc=0 of the Brillouin zone they sum up to 
the value: 4:{Vpda - 2/%/3Fpd7r)- 



FIG. 2: Mnoa— Mni pair in the GaAs structure. 

In contrast, the ion in a tetrahedral interstitial posi- 
tion, e.g., as in Fig. 0, has 6 anion neighbors on 
[0, 0, ± 1], [0, ± 1, 0] and [±1, 0, 0] directions, at the 
distances a/2 (see Fig. |). In this case all not equal zero 
inter-atomic matrix elements are proportional to the ap- 
propriate sm{aki/2) (where ki, i — x,y,z, are the com- 
ponents of the wave vector k) and they vanish at the 
centre of the Brillouin zone. Thus, the Mnj d-orbitals 
do not hybridize with the p-states of the holes at the 
top of the valence band (i.e. for the tetrahedral intersti- 
tials the kinetic exchange constant Pi^— 0) and they do 
not contribute to the hole-induced ferromagnetism. This 
means that the formation of Mn interstitials decreases 
not only the hole concentration but also the number of 
Mn ions participating in the Zener-type ferromagnetism. 
Still, these effects do not explain why the removal of in- 
terstitials leads to the increase of magnetization and to 
the higher Tg. than expected from the rise of the hole 
concentrationo. j-. 

As pointed out already by Yu et alE the electrostatic 
attraction between positively charged Mni donors and 
negative MuQa acceptors stabilizes the otherwise highly 
mobile Mnj in the interstitial sites adjacent to MuQa, 
forming a Mnca-Mni pair, as shown in Fig. |^. One no- 
tices that despite the fact that for the interstitials the 
p-d kinetic exchange and, consequently, the hybridiza- 
tion mediated spin interactions with holes in the vicinity 
of the top of the valence band vanish, the ionic spins in 
the pair can be coupled by superexchange mechanism. In 
the latter process the spins of the two ions. Si and 5*2, 
are correlated due to the spin-dependent p-d exchange in- 
teraction between each of the ions and the valence band 
electrons in the entire Brillouin zone. The superexchange 
Hamiltonian: 

H superexchange — 

-2J(i?i2)5i • ^2 (1) 

can be obtained by a proper selection of spin-dependent 
terms in the matrix of the fourth order perturbation with 
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respect to the hybridization for a system of two ions in 
the crystal: 



^ {l\H,,\l"){l"\H,,\n{V\H,,\l){l 
1,1', I" 



(2) 

Using the virtual transition picture, one can say that 
the superexchange is a result of four virtual transitions 
of an electron - from the band onto the d-shell of the 
ion and fmm the ionic d-shell to the band, in different 
sequencescJ. The quantitative determination of the su- 
perexchange constant J requires the knowledge of the en- 



ergies of these virtual transitions, which are represented 
by the energy differences between the intermediate and 
initial states of the system of two ions and the completely 
filled valence bands, in the denominator of Equation (2). 
Of primary importance it is, however, to determine the 
sign of the superexchange interaction for the Mnca-Mni 
pair. In the following, we calculate the exchange con- 
stant J within a simplified model, in which we neglect 
the dispersion of the valence bands but we account for 
the wave- vector dependence of the hybridization matrix 
elements. The resulting formula for the exchange con- 
stant J reads: 



E„ 



EliEa,+E,,) EliEa,+E,,) 



(3) 



V V* r {'^W % (2)y* r (XW r (1) 



vi ,1^2 ,ki ,k2 ,m ,n 



In Equation (3) the summation runs over the valence 
band indices vi and i>2, the wave- vectors ki, k2 from the 
entire Brillouin zone, and over the Mn d-orbitals m, n. 
The energies Ea^ and Ed^ {i = 1,2) are the transfer en- 
ergies for the electron from the valence band onto the 
ion i (" acceptor" ) and from the ion i to the valence band 
("donor"), respectively. It should be noted that these en- 
ergies for the interstitial Mn ion are completely unknown. 
Still, since all these energies as well as the sum, which 
we calculated numerically, are positive, we can conclude 
that the Mnca-Mni pair is antiferromagnetically coupled. 
Thus, Mn ions when in tetrahedral interstitial positions 
not only do not contribute to the hole-induced ferromag- 
netism but they also make some of the substitutional Mn 
ions magnetically inactive by forming with them close 
pairs, in which the spins of the ions are antiferromagnet- 
ically coupled by the superexchange mechanism. This 
explains the experimental observations that the removal 
of Mni ions by low-temperature annealing leads not only 
to an increase of the hole concentration, but also to a 
significant increase of the magnetization. 

To estimate the strength of this coupling we compare 
J with the superexchange constant J' for a Mnca-Mnca 
closest pair, obtained within the same simple model. 
Using the same transition energies for both Mnca and 
Mni ions, we obtain J/ J' « 1.6. This is not surpris- 
ing in view of the small distance between the interstitial 
and the nearest substitutional Mn ions and the larger 
number of anion neighbors for Mnj. With a reasonable 
value of 3 eV for the Mnca charge transfer energies, with 
the values of Harrison parameters Vpda=^-^ eVE3, and 
VpdTT = —\Vpd(j, (which for tlmMni we scale according 
to the Harrison's prescription)^ the absolute values of J 



and J' constants are by far not negligible: J w 71 K and 
J' w 43 K. 

The role of interstitial Mn ions occurred to be even 
more pronounced in Gai-j^-yMn^jBeyAs samples, grown 
at Notre Dame with the hope to increase the hole 
concentr|ation, and hence Tc, by introducing another 
acceptorcS. Instead, it turned out that adding Be to 
Gai_a;Mn2;As inpieases the concentration of Mnj at the 
expense of MnGaQE^. At the same time, although the hole 
concentration (joes not change significantly, the Tc drops 
dramaticallyMES, in agreement with the presented above 
result that the Mnj do not participate in the hole-induced 
ferromagnetism. Recently performed angular scans seem 
to suggest, however, that the Beca acceptor stabilizes 
the Mn interstitial donor not in— the tetrahedral but in 
the hexagonal, ( ||| ) positionc2l. In this site the Mni 
has three anion nearest neighbors, as shown in Fig. H, on 
the [-1, -3, -1], [3, 1, -1] and [-1, 1, 3] directions, at the 
distance a-s/TT/S. In such case one does not expect the ki- 
netic p-d exchange constant (3 to be equal to zero - it can 
be rather expected that the hybridization for the ion in 
this site should be stronger, due to the smaller distance to 
the anions, what increases the Harrison parameters Vpdu 
and Vpdir)- Surprisingly enough, the most of the inter- 
atomic matrix elements in the hybridization constant V 
for the hexagonal interstitial Mn mutually cancel at the 
centre of the Brillouin zone. This leads to a consider- 
ably smaller than for the substitutional Mn ion value of 
the kinetic exchange constant /3/^, i.e., ~ 5. As the 

Curie temperature in the Zener model depends on 0^ , we 
conclude that the contribution to the hole-induced ferro- 
magnetism from the Mn ions occupying the hexagonal 
interstitial sites is as well very much suppressed. 
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FIG. 3: The six (three cations and three anions) nearest 
neighbors and the next four cations and four anions for an 
ion in the hexagonal interstitial position in the zinc blende 
lattice. 



In conclusion, we have shown that not only the com- 
pensating properties of the interstitial magnetic ions 
impose a limit to the Curie temperature in the fer- 
romagnetic Gai_2:Mnj:As and Gai_2:_j^Mna;Bej,As sam- 
ples. Also their magnetic properties in both (tetrahe- 
dral and hexagonal) interstitial sites, i.e., the negligible 
kinetic exchange constant and strong antiferromagnetic 
superexchange with the adjacent substitutional Mn ion, 
act towards diminishing the transition temperature. 
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